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Abstract 

For atomic thin layer insulating materials we provide an exact analytic form of the two- 
dimensional screened potential. In contrast to three-dimensional systems where the macroscopic 
screening can be described by a static dielectric constant in 2D systems the macroscopic screening 
is non local (q-dependent) showing a logarithmic divergence for small distances and reaching the 
unscreened Coulomb potential for large distances. The cross-over of these two regimes is dictated 
by 2D layer polarizability that can be easily computed by standard first-principles techniques. The 
present results have strong implications for describing gap-impurity levels and also exciton bind- 
ing energies. The simple model derived here captures the main physical effects and reproduces 
well, for the case of graphane, the full many-body GW plus Bethe-Salpeter calculations. As an 
additional outcome we show that the impurity hole-doping in graphane leads to strongly localized 
states, what hampers applications in electronic devices. In spite of the inefficient and nonlocal 
two-dimensional macroscopic screening we demonstrate that a simple k • p approach is capable to 
describe the electronic and transport properties of confined 2D systems. 

PACS numbers: 73.22-f, 78.67.-n, 71.35. Cc 
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INTRODUCTION 



The study of two-dimensional (2D) electronic systems is of great fundamental significance 
in physics. Atomic thin layers allow to address the role Coulomb interactions in confined 
geometries [lj]. In this context the synthesis of graphenej^ has triggered a huge amount of 
work in the understanding and controlling the properties of this system. In fact due to its 
unique electronic properties and the low dimensionality, graphene is considered as one of 
the most promising materials for future carbon-based electronics. Nevertheless, the peculiar 
gapless ultra- relativistic energy spectrum of graphene 0, makes the creation of carbon 
nanodevices based on p-n junctions highly nontrivial. Therefore, transforming graphene 
into a semiconductor with a conventional electron spectrum keeping its two dimensionality 
introduced a challenge that has been a major line of research in the last years. Recently 
an important step towards graphene electronics has been made with the synthesis of a 

functionalized 

Graphane is a wide band-gap dielectric [ 



fu.lv hydrogenated graphs™ named g r apha„e H as we.l as other chemical 
graphene-like structures [6|, [7J 



9| and therefore 



it may become an important part of nanoelectronic devices as it opens a way to create 



2D p-n junctions 
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llllj . Similarly, single-atomic layers containing hybridized domains of 
121 ] have been synthesized and follow a completely different and novel 



graphene and h-BN 

(electronic) phenomenology as compared to high-purity or damaged graphene. Actually, the 
search for low-dimensional semiconductors is not only focus on graphene and its derivative 
compounds but it is also moving towards other layered systems as for example M0S2, WS2, 
MoSe2, MoTe2 and BN which can be efficiently dispersed in common solvents and can be 
deposited as individual flakes or formed into films 13|. In fact, M0S2 monolayer has been 



now synthesized 



14 



band gap of 1.8 eV 



3- 



Contrarily to the bulk M0S2, it is a direct gap semiconductor with a 



14J. Thus, it is important 



and could be used as single-layer transistor 
at this point to make a deep analysis of many body effects and in particular of the nature 
and functional form of the screening in general low dimensional systems (in particular two- 
dimensional semiconductors and insulators). In fact, screening effects play a fundamental 
role in determining the electron dynamics, the exciton binding energy and the effective 
electron-electron and electron-phonon interactions in the superconducting state. Moreover 
the screening dictate the optical and transport properties of 2D devices so that knowing its 
behavior in low dimensional systems is fundamental also for practical applications. 



systems has been 



18| . to the best of 



Although in the past the problem of the screening in low dimensiona 
extensively investigated in 2D metals [l7] and semiconducting thin films 
our knowledge it has never been dealt in a strict 2D dielectric that will be the focus of this 
letter. In the present work we provide a strict 2D derivation of the macroscopic screening 
derived by Keldysh as a limiting case of a thin filml8|. We demonstrate that, contrarily 
to what happens in 3D systems where the macroscopic screening is mapped in a dielectric 
constant, in 2D systems the macroscopic screening is non local so that in the Fourier space 
it is described by a q dependent macroscopic dielectric function. 

Among the various 2D dielectrics, graphane is not only a promising material for nano 
devices application, but also is very interesting by itself. In fact, theoretical works based on 
first principles calculations predict localized spin states at hydrogen vacancies 19[ , demon- 



transfer excitons 



201 ] and indicate that 



211 ] . Therefore in the present work we 



strate the existence of unusual strongly bound charge- 
doped graphane is probably a high T c superconductor 
take graphane as the test system to address in detail the influence of the 2D screening of 
the Coulomb potential on the excitonic states and impurity levels. Our findings are general 
in scope and can be applied to any other 2D insulator as the ones described above. 

The present work is organized as follows. First, we derive in a simple electrostatic model 
the exact two-dimensional screened potential and compare it with its three-dimensional 
counterpart. We also provide a very simple and pictorial understanding of that potential 
in terms of the potential created by a one-dimensional charge distribution which length 
is determined by the two-dimensional layer polarizability. The effect of the specific 2D 
macroscopic screening is illustrated addressing the electronic and optical properties of perfect 
graphane as well as the electronic levels introduced by hydrogen vacancies in the layer (hole 
doping). We close the paper with some brief conclusions and perspectives. 



DIELECTRIC SCREENING IN 2D INSULATORS 

To determine a long wavelength static dielectric response of a general 2D insulator we 
consider a dielectric sheet of zero thickness at z = embedded into vacuum, and subject to 
an external potential </>ea;t( r )- For definiteness we assume that <f) e xt{^) is produced by a point 
charge placed at the origin n ext (r) = e5(r). The total electrostatic potential <p produced by 
the external source is related to the total charge density n by Poisson's equation: 
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V>(r) = -47m(r) (1) 

where n = n ext + Bmd is the sum of the external charge density and the induced charge 
density. The induced charge density is confined on the plane z = and, in the long wave 
length limit, is related to the 2D macroscopic polarization ~P 2 d {jimd — — V-P2d)j which, in 
turn is proportional to the in-plane component of the total electric field. Introducing the 2D 
polarizability a 2 D of the dielectric sheet, so that P 2 £>(p) = — «2D V p 0(p, z = 0) we obtain 
an expression of the induced charge density in terms of the macroscopic potential evaluated 
at a point r = (p, z — 0): 

rw(r) = S(z)a 2D V 2 p (/){p, z = 0) (2) 

With this result for the induced charge density the Poisson equation for the potential of the 
external point charge takes the form: 

V 2 0(r) = -47re5(r) - 4tt« 2D V 2 0(p, z = 0)6 (z) (3) 

while its Fourier transform can be written as: 

/elk 
Z7T 



<^2D(q) 

where q is the in-plane component of the wave vector, and the /c^-integral in the right hand 
side defines the Fourier component </>2d(q) of the 2D macroscopic potential. By solving 
Eq. (@J we obtain the following result for <^2r>(q) 

02d(q) = | ft; 2 ^ e rT\ ( 5 ) 

|q|(l + 27ra 2 D|q|) 

which defines the 2D macroscopic screening of a point charge. As can be seen, for a 2D 
insulator the macroscopic dielectric screening is no more described by a simple dielectric 
constant which renormalizes the electronic charge as in 3D systems. In 2D systems a formally 
defined dielectric function is intrinsically q-dependent: 

e(q) = 1 + 27ra 2j D|q| (6) 

By the inverse Fourier transform of e0 2 £>(q), we can determine the effective potential V e ff(p) 
which is felt by an electron living in the 2D dielectric in presence of a point charge: 



Veff(p) 



4a 



2D 



Y 



(7) 



where H and Y are the Struve function and the second kind Bessel function respectively and 
r = 2ira2D- From the known asymptotic properties of the Struve and Bessel functions j^J 
we determine the following asymptotic behavior of V e ff(p): 



HmV.„(,)~± (8) 



lim V e ff(p) ~ 



(9) 



In ( — — ) + (7-ln2)e £ 



(10) 



where 7 ~ 0.5772 is the Euler's constant. By the simplest possible matching of the two 
asymptotic behaviours we can construct an approximated expression for V e ff(p) in terms of 
elementary functions: 

,P + r o, 

which gives an accurate description of the effective interaction also at intermediate values 
of p/r as can be inferred from FigJT] 

The above results clearly show that, in contrast to the 3D case, the screening in 2D 
dielectrics introduces a new length scale r$ which is determined by the the polarizability 
«2D of the dielectric layer. When p is larger than r the effective potential behaves like the 
3D unscreened Coulomb potential while for p — > it diverges logarithmically, i. e., it goes 
like the Coulomb potential in two spatial dimensions (the potential of a charged string). 
Obviously the logarithmic divergence weakens when a2D increases, which means that the 
screening is more efficient in highly polarizable systems. 

To better understand the difference of the screening in 3D and 2D insulators we consider 
a point charge surrounded by a 3D and 2D dielectric medium, respectively (see Figj2] (a) 
and (b)). The total electric field (E) at a distance r from the point charge will be the sum 
of the external field produced by the point charge (E ext (r) = ^f) and the induced field 
(Ej n d(r) = — 47rP(r)). In 3D dielectrics the latter is equivalent to the electric field produced 
by a uniform charge distribution on a sphere of radius r centered on the point charge (see 
Figj2] (a)). This charge distribution produces a field of the same functional form as that of 
the external point charge itself, E inrf (r) ~ E e:r4 (r), which means that the screening is given 



by a simple multiplicative renormalization. In the 2D case the situation is quite different. 
As can be inferred from Figj2] (b) , since the system is polarizable only on the plane, Ej n d is 
equivalent to the electric field produced by a uniform charge distribution on a circle of radius 
r. As a consequence it will be a function of r and 9 with a functional form substantially 
different from E exi (r). This results in a non-local macroscopic screening. 

A simple and pictorial understanding of the 2D effective potential V e ff(p) can be obtained 
by rewriting Eq. ([7]) in a different form. Starting from eqj5] we replace the factor (1 + 
27ra2Z)|q|)~ 1 by its integral representation and rewrite 02i)(q) as follows 

i*i 

02i)(q) = / dz—e |q||z| - . (11) 

J-oo |q| 47To;2D 

Performing the Fourier transform we get the following expression for the effective interaction 
between an electron and an external point charge 

_ |z| 

v -"^ = L dz 7FT7^ (12) 

Obviously this equation is absolutely equivalent to Eq. (jJJ), but it is much more clear phys- 
ically. Indeed, Eq. ( fT2|) represents the potential (in the plane z = 0) produced by a one- 
dimensional charge distribution of the form 

_M 

Q(r) = e5(p) e —^ (13) 
zr o 

Noting that / Q(r)d 3 r = e, we conclude that in the presence of the dielectric plane the point 
charge produces a field as it would be effectively smeared out into a ID string with the 
charge distribution of Eq. ( [TBI . This behaviour should be contrasted to the multiplicative 
renormalization of the charge in 3D dielectrics. 

Thus, the effect of the 2D dielectric screening can be visualized as follows: two electrons 
living in a 2D dielectric plane interact as two thin charged rods of the length ~ 2tq and 
the line charge density Q(r) ~ e/2r . The length of the rod sets the characteristic scale 
of the potential. From large distances p » the rod is seen as a point charge with 
the potential given by the classical 3D Coulomb law Eq. (JED. Hence at large distance, the 
induced polarization is completely inefficient in screening the external field. In the opposite 
limit p « r the rod looks like an infinite wire with the line charge density ^ so that the 
effective potential reduces to the classical 2D Coulomb potential of Eq. Q. Thus at small 
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distance the effect of the induced polarization becomes dominant - the 1/r singularity is 
replaced by a weaker logarithmic dependence. 

It should be noted that our results for the 2D dielectric screening are very closely related 



to the results obtained by Keldysh in ref. 



for the interaction potential of two point 



charges in a dielectric slab of the thickness d and characterized by a static bulk dielectric 

and references therein). In fact, our Eq. (J7J) can be recovered in 



constant e (see also Ref. 
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the limit p> d and e> 1 [18]. The ID distribution of the effective charge Eq. ([13]) can be 
also viewed as a limiting form of the discrete image charges used to construct the solution 
of the electrostatic problem for a finite dielectric slab 23j. The important novel outcome of 
our derivation is that the form of the effective screened potential of Eq. (J7J) is valid even for 
a microscopically 2D, atomically thin dielectrics for which the notion of the bulk dielectric 
constant makes no sense. 

The only parameter entering the screened potential of Eqs. (0) and (J7J) is the polarizability 
«2D of the 2D dielectric. Let us show how it can be extracted from the standard ab-initio 
supercell calculations where 2D systems are simulated using a periodic stack of layers with 
sufficiently large inter-layer distance L. For this auxiliary 3D layered system we can get the 
3D macroscopic polarization P 3 e> = a^^E, where a^o and E are the 3D polarizability end 
the total electric field respectively. The macroscopic 3D polarization can be calculated as 
an average over N layers in the periodic stack of the microscopic 3D polarization P^ c (z): 

P ^ = ^ZJ dzP 3D C (z) (14) 

where, in the definition of P^ c {z) we have already performed the one-layer average. Hence 
P^ c (z) can be expressed in terms of the macroscopic 2D polarization P2D = <^2dEi oc as 
follows: 

N 

Pyoiz) = E P 2dS(z - nL) (15) 

n=0 

with E[ oc being the local field acting on a single layer. Inserting Eq. (TI5]) in Eq. (1141) and 

taking L sufficiently large so that E\ oc rs E we obtain an expression of a 2 D in terms of a^D 

(as a generalized Clausius-Mossotti expression for 2D systems): 

e - 1 . . 

a 2 D = La 3D = L— — (16) 

47T 

where e is the static dielectric constant of the 3D layered system. The value of e entering 
Eq. ([16]) can be evaluated directly from the first principles calculation of the dielectric 



function €gg'(<1,u) as: 

e = lim - — — ^ . (17) 

q^o [e-ifau = 0)] g =g'=o 

When L goes to infinity, e of the layered system approaches the vacuum dielectric constant, 
e = 1 + 0(1/L). Therefore the L — > oo limit of the right hand side in Eq. (fT6l) yields a finite 
value that is equal to the 2D polarizability. In practice one performs calculations for several 
sufficiently large L to ensure the convergence of a>2D- 



EXCITONIC AND IMPURITY STATES IN GRAPHANE 



We apply the results of the previous section to the description of excitonic and impurity 



jeen addressed recently using 
In this section we show 



states in graphane. The problem of excitons in graphane has 
a fully ab initio many-body self-energy GW-BSE approach 
that the effective screened potential of Eq. ([7]) combined with the k • p description of the 
electronic and hole states leads to a very simple and accurate description of strongly bound 
electron-hole and hole-impurity states as obtained from the GW-BSE calculations. 

Graphane is a representative of wide band gap 2D dielectrics, which is obtained from 
the ideal graphene by depositing hydrogen atoms on both sides of graphene plane. The 



resulting electronic structure is dictated by sp 3 hybridization of t 



causes the opening of a wide band gap (of 5.4 eV at the T point) js|. States at the top of the 
valence band belong to E 9 2D irreducible representation of the graphane point group T)3d, 
while the bottom of the conduction band belong to the A 2u ID irreducible representation. 
As shown in Ref. [2^ transitions from the top of the valence band to the bottom of the 
conduction band are allowed in the dipole approximation and result in strong excitonic 
effects in the absorption spectra. In particular, the corresponding electron-hole pairs give 
rise to two nearly degenerate excitons with binding energy of about 1.6 eV. This large 
binding energy (one order of magnitude larger then in typical semiconductors) seems to be 
surprising since both valence and conduction bands form almost perfect parabolas in a wide 
energy-momentum range around the T point and therefore the excitonic states are expected 
to be well described in terms of the effective mass approximation, in spite of their small 
radius. We will demonstrate explicitly that the effective mass approximation does indeed 
works perfectly, and that the unusually large binding energy is completely explained by the 
weak and nonlocal 2D screening discussed in the previous section. 
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re carbon orbitals, which 



Angular 


Quantum 


Exciton 


Impurity 


momentum (I) 


number (n) 


energy (eV) 


level (eV) 


I = ±1 


n = 1 


-1.77 


-2.12 


n = 2 


-0.67 


-0.90 




n — 1 


E + = -1.13 


E + = -1.47 


I = 


= -0.78 


E_ = -0.92 


n = 2 


E + = -0.52 


E + = -0.74 




= -0.34 


E_ = -0.45 


I = ±2 


n = 1 


-0.92 


-1.17 


n = 2 


-0.43 


-0.58 



TABLE I. Exciton binding energy and impurity levels in the effective mass approximation for some 
selected values of the quantum numbers n and I. 
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Let us start with the k-p effective mass approximation for the electronic states in graphane 

2/)i,. 



24]. The Hamiltonian for the conduction band is trivially given by H c (p) = ^2rn V ' wnne 
for the valence band Hamiltonian H v (p) we adopt the representation obtained in Ref. 241 ]: 

(18) 



H v (p) = l -alp 2 + ^[a + pl + a_tf_ 



where p = — iW is the in-plane momentum operator, p± = p x ±ip y , I is the identity matrix, 
and a± — g x ± io y with aj being the Pauli matrices, and a = 2.62/mo, /3 = 0.98/m and 
m e = 0.83m are the band parameters expressed in terms of the bare electronic mass m 
and obtained from the ab-initio band structures 201 ] . 

The excitonic Hamiltonian for the zero momentum excitons can then be constructed in 
standard way 25]: 

H ex = H c (p) + H v (p) - V eff (p) (19) 

where V e ff(p) is the effective 2D screened electron-hole interaction given by Eq. (J7J). Ex- 
plicitly the final effective mass equation for the relative motion of the electron and the hole 
takes the following form 



1 - 1 

^/Jp 2 + ^72[<3-+P+ + fr-pl] - V eff {p) ) $(p) = E${p) 



(20) 



where 7i = « + and 72 = P- 



To classify the eigenstates of Eq. f[2"0"j) we note that the Hamiltonian H ex commutes with 
an operator L z that is defined as follows 



IP x P) 



a. 



(xp v -ypx)- cj z 



(21) 



Obviously, the operator L z corresponds to the ^-component of the total angular momen- 
tum, with the second term in Eq. (T2~T]) being related to the orbital momentum of the local 
currents inside the unit cell of graphane 24|. Since [H ex , L z ] = the excitonic states can be 
classified by the eigenstates of the total angular momentum operator. In other words, the 
eigenf unctions of Eq. (12"U|) can be written in terms of radial wave functions (Zi(p), Xi(p)) 
ordered by the integer quantum number I defining the eigenvalue of L z : 



( Ad 



\ 



i6 Xi{p) 



(22) 



Inserting the expression of Eq. 
part of the envelop wave functions: 



into Eq. (|2"U|) we obtain the equation for the radial 



7i 
' 2 



p p p 



7i 
' 2 



dl + -d p 



p , 
[i-if 



-K//(p)J^(p)- 

Xi(p) = E&ip) 
Zi{p) = E lX i{p) 



(23) 



(24) 



Thus, each excitonic state is completely defined by the quantum number / and the positive 
integer n denoting the discrete eigenvalues of Eqs. f )23j) and (T24"j) for given I. The correspond- 
ing microscopic wave function of the exciton for a fixed position of the hole can be written 
as follows: 



^(r,r,) = Z i (p)^(r)^ 1) (r,) + 



(25) 



where ip e is the electron Bloch wave function and %l) h ' the hole Bloch wave functions related 
to the two fold degenerate valence bands. 
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Analyzing the structure of Eqs. (|23|) and (|24|) . we observe that for all I ^ the system 
of differential equations is invariant under the transformation / — > —I, Z — > x, \ ~ * 
Therefore all excitonic states with / ^ are double degenerate with Ei = E-i, which is a 
clear consequence of the time-reversal invariance of the Hamiltonian. We also note that only 
excitons corresponding to I = ±1 are dipole active. 

The only non-degenerate state corresponds to a dark exciton with zero angular momen- 
tum, / = 0. Interestingly, for / = the diagonal and off-diagonal operators in the system 
of Eqs. (j23|) . ([21]) are equal to each other. As a results the problem reduces to completely 
decoupled equations for the "symmetric" and "antisymmetric" states 



±72) U + U P ~ ± J - V eff (p) 



MP) = E ± MP) (26) 
and the excitonic spinor wave function for I = takes the form 



1 / e id cf>± 



Co = 4= ,„ (27) 



y/2 \ e- id ± 



To practically solve Eqs. ( 123]) and (I24j) we expanded the radial part of the envelop wave 



function on the 2D hydrogen eigenf unctions f[26|. so that Zi(p) = Y, n a nU n ,i+i(p) and 
Xi(p) = J2nb n u n ,i-i(p)- This complete orthonormal basis set assures the correct asymptotic 
behavior of the eigenfunctions of the excitonic Hamiltonian. 

Our results for a selected set of lowest energy states are summarized in table [B As we can 
see, the ground state of the excitonic Hamiltonian corresponds to I — ±1. This state is two- 
fold degenerate and optically active with the binding energy Ei = ±i = 1.77 eV, which is in a 
perfect agreement with the values obtained by solution of the Bethe-Salpeter equation j^. 
The corresponding excitonic wave functions (Figj3]), calculated using Eq. (1251) . shows that 
both excitons are strongly localized with an average radius of about 11.5 a.u. As can be 
inferred from Figj3] these excitations give rise to a charge transfer from the carbon plane 
towards the hydrogen plane. The first excited state correspond to zero angular momentum 
I = 0. This dark exciton is also found from the solution of the Bethe Salpeter equation, 
which is additional confirmation of the present simple theory. 

Therefore our results demonstrate that, despite the large binding energy, excitons in 
graphane are indeed described in terms of the effective mass approximation, provided the 
correct form of the effective electron-hole interaction is used (as derived in the present work, 
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Eqs. flTJ) and (fT2|) . We clearly see that the unusual, large binding energy is related to a weak 
and nonlocal 2D dielectric screening which is completely inefficient at large distances. The 
small overestimation of the exciton binding energy respect to the ab-initio value may be 
ascribed to the lack in our approach of short range contributions to the induced polarization 
and exchange electron-hole interaction. All these effects can only reduce the exciton binding 
energy. As a matter of fact, the effects of short range corrections are small and, if necessarily, 
can be easily included perturbatively. 

Using the same formalism we can now look at the effect of the 2D screening on impurity 
states. First, we focus on the acceptor states in the hole doped graphane, as it is expected 
to be the most natural way to dope this system. Indeed naively one may assume that extra 
holes are easily introduced by dehydrogenation. In this context it is worth noting that hole 



doped graphane has been predicted to be a high T c superconductor [211 ] . 

Similarly to the excitonic case, the parabolicity of the valence bands in a wide energy 
range suggests that holes in the presence of hydrogen vacancies can be well described in terms 
of the effective mass approximation. Therefore acceptor impurity levels can be obtained by 
solving Eqs. (I23p . ff24|) with 71 = a. The corresponding results are presented in table [B 
The ground impurity state corresponds to n = 1 and I = ±1 and is characterized by a 
binding energy of about 2.12 eV in good agreement with the ab-initio value (1.86 eV). This 
quantity represents the position of the impurity level with respect to the top of the valence 
band. Therefore, for impurity levels the 2D nonlocal screening results in unusual large 
binding energy that exceeds by two or three orders of magnitude the corresponding values 
for typical semiconductors. Comparing the values of the binding energy with the gap energy 
(5.4 eV) we find that the impurity level is close to the center of the graphane gap. 

Finally, when the dopant is a donor, the electron in the conduction band is described by a 
simple 2D hydrogen like Schrodinger equation with K//(p) of Eq. ([7]) replacing the Coulomb 
potential. In this case for the lowest bound state (corresponding to n = 1 and / = 0) we get 
a binding energy of about 3.15 eV. Therefore for electron doped graphane the specific 2D 
screening of the impurity potential also causes the formation of mid-gap impurity levels. 

The above results lead us to an unfortunate but important conclusion. The standard 
for 3D semiconductors impurity doping, both donor and acceptor, will probably not work 
for graphane and most likely for other atomically thin dielectric. In particular a slightly 
dehydrogenated graphane cannot be considered as a semiconductor with extra highly mobile 
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holes in the valence band. All holes will be strongly localized on the hydrogen vacancies 
with the radius of the bound state of the order of the interatomic distance. The reason for 
this behavior is a very weak and inefficient screening in 2D dielectric materials. 

CONCLUSIONS 

In conclusion, we derived an expression of the macroscopic screening in 2D dielectrics 
showing that, contrarily to what happens in 3D systems where the macroscopic screening 
is mapped in a dielectric constant, in 2D systems the macroscopic screening is non local. 
As a result the effective potential produced by an external point charge surrounded by a 
2D dielectric has a functional form which is substantially different from the bare Coulomb 
potential. It presents a logarithmic divergence for p — > and reduces to the unscreened 
Coulomb potential at large distances. The 2D polarizability a^n determines the character- 
istic length scale tq at which the two asymptotic forms are matched. This behavior strongly 
modify the optical and transport properties of 2D systems. In particular we show that hole 
impurity doping leads to strongly bound localized states with low mobility. Moreover, spite 
of the inefficient and two-dimensional macroscopic q-dependent screening the simple k • p 
approach works very well to describe the electronic properties up to very high energy, and 
very short spatial scales. Our results imply that the k ■ p theory supplemented with a proper 
macroscopic treatment of the 2D screening forms a solid basis for a quantitative descrip- 
tion of various, both equilibrium and nonequilibrium, in particular, transport properties of 
nanostructured 2D systems. 
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FIG. 1. Comparison between the true effective potential V e ff(p) from Eqf7]and its approximated 
form Vejf(p) described by EqfTOl 
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FIG. 2. Schematic representation of the effect of the macroscopic polarization induced by a positive 
point charge on the z=0 plane in 3D (a) and 2D (b) dielectrics. 
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FIG. 3. 3D-Shape of the low energy excitonic wave functions for a fixed position of the hole (marked 
as a green circle) as obtained from Eq. (|25p . Note that the shape of the excitonic wave functions 



is in perfect agreement with that obtained by the full solution of the BS equation in Ref. 
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